Fracture Roughness and Correlation Length in the Central Force Model 
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We measure the roughness exponent and the correlation length exponent of a stress-weighted 
percolation process in the central force model in 2D. The roughness exponent is found to be £ = 
0.75 ±0.03 and the correlation length exponent is found to be v = 1.7 ±0.3. This result supports a 
conjecture that the fracture roughness for large scales is controlled by a stress weighted percolation 
process, and the fracture roughness can by calculated from the correlation length exponent by 
£ = 2z//(l + 2u). We also compare global and local measurements of the fracture roughness and do 
not find sign of anomalous scaling in the central force model. 

PACS numbers: 62.20.Mk 
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In the early eighties, it was observed that brittle frac- 
ture surfaces show self affinity. That is, such a surface 
will be statistically invariant if the in-plane length scales 
are changed by a factor A, and the out-of plane length 
scale is changed by a factor AS where C is the rough- 
ness or Hurst exponent In the early nineties, it was 
proposed that the roughness exponent is universal Q. 
This initiated a large effort to further investigate this 
henomenon An early review may be found in Ref. 



For large scales the universal value of the roughness 
exponent is believed to be about 0.80 for brittle fractures. 
There is also evidence that there may be a different value, 
0.5, for small scales, with a well-defined crossover length 
scale separating these two regimes 0- 

Even though the evidence for one or more universal 
roughness exponents is mounting, the mechanisms that 
may be responsible are still not known. There have been 
some proposals for mechanisms, see 0, 0, H, @- It is 
the aim of this Letter to test the basic mechanism pro- 
posed by Hansen and Schmittbuhl connecting the 
roughness exponent to a correlation length exponent as- 
sociated with the underlying fracture process. The basic 
idea is that the breakdown process when the material 
is very disordered is essentially a correlated percolation- 
like process — percolation like in the sense that the ma- 
chinery for describing it is the same as in ordinary per- 
colation. For smaller disorders, the breakdown process 
localizes, and by combining the localization length with 
the percolation-like description, a relation may be set up 
between the roughness exponent £ and the correlation 
length exponent v of the correlated percolation process, 



c = 



2;/ 



\ + 2v 



(1) 



This relation has been tested on the two and three- 
dimensional fuse model [Tfl ITlj . The fuse model con- 
sists of a regular network of electrical fuses, each having 
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a burn-out threshold drawn from some spatially uncor- 
related statistical distribution [l2L fl3| . The model has 
been used quite extensively to study fracture roughness 
H III Elll El In this Letter, we test Eq. JTJ in the 
central-force breakdown model, first studied in Ref. [Isj . 
This is a model that is much closer to the fracture prob- 
lem in that its response is elastic — in the continuum 
limit it maps onto ordinary Lame elasticity. It consists 
of a regular lattice whose bonds are elastic springs free 
to rotate around the nodes they are connected to. The 
force on a spring between nodes r% and r*j given as |Tflj 



h 



( r 3 - n)) 



(2) 



where cry is the spring constant for the spring and n is 
the axis vector of the spring. Another elastic model that 
has been studied is the beam model [2fl l2ll|. 

Measuring both £ and v independently in the two- 
dimensional central-force breakdown model, we will be 
able to test Eq. for this case. For pure central-force 
rigidity percolation in two dimensions, Moukarzcl and 
Duxbury found v = 1.16 using the pebble game al- 
gorithm. As is the case for the fuse model 0,H3, there is 
no reason whatsoever that the v of the breakdown process 
should be equal to the percolation correlation exponent. 
In the fuse model, one finds v = 1.56 in two dimensions 
to be compared to v — 4/3 in two-dimensional percola- 
tion, and v = 0.83 ± 0.04 in the three-dimensional fuse 
model, to be compared to v = 0.88 for three-dimensional 
percolation. As we shall see, we find v = 1.7 ±0.3 for the 
two-dimensional central-force breakdown model. 

We simulated fractures in 2D trigonal central force lat- 
tice with periodic boundaries in the x-direction and ap- 
plied tension in the y-direction. Each spring in the lattice 
was assigned a threshold i; at random from a power law 
distribution 7f(i) oc which gives ti — rf where ri is 

a random number between and 1 and D = (3 . Values 
of D close to zero gives narrow disorders. The disorder 
becomes broader when D increases. We used D = 0.7 
for the narrow disorder and D = 20 for the broad disor- 
der. For each realisation, a strain of unity is applied in 
the vertical direction to the boundary lattice points and 
the central force equations are solved iteratively with the 
conjugate gradient method [24j, The lattice is periodic 
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FIG. 1: The global roughness exponent is found by consid- 
ering the scaling of the average fracture width of a fracture 
profile y(x), using Eq. where the averages are done over 
the whole system. The solid line is £ = 0.75. We used from 
5499 samples for L = 16 to 237 samples for L = 256. 

in the horizontal direction. The spring with max^ (/;/£;) 
is identified and the strain is increased (decreased) until 
fi /U — 1 and the spring constant for this spring is set to 
zero. The central force equations are solved again with 
this new configuration and springs are removed until the 
elasticity module of the system is zero. 

When a sample is fractured, elasticity module equal 
zero, the fracture surface is created in the dual lattice 
which gives 2 x L x measuring points for the fracture sur- 
face for the triangular lattice. Before the fracture width is 
calculated, overhangs are removed using a solid-on-solid 
approximation. 

In order to address the question of whether there 
is a possibility for anomalous scaling of the roughness 
[25L I2r1 |. we distinguish between measuring the rough- 
ness exponents using global methods — that is, when the 
roughness is measured along the entire length of the sys- 
tem, and the system size is changed, and local methods, 
where the size of a window is changed whilst the system 
size is kept fixed. When global and local methods give 
different roughness exponents, the system is anomalously 
rough. 

We calculate the global fracture roughness exponent 
by finding how the fracture width scales with system size 

W(L) = ((y(xf) L -{y(x)))iy/ 2 <xLS , (3) 

and the local fracture roughness exponent with the scal- 
ing of the local width with window size 

w(l) = ((y(x) 2 ) l - (y(x)W 2 oc , (4) 

and the average wavelet coefficients (AWC) method [27| 

W[y](a) oc a Ci ° c+1/2 . (5) 

The simulations for measuring the fracture roughness 
was done with a narrow disorder of D = 0.70. For the 
global exponent we measured £ = 0.75 ± 0.03. 
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FIG. 2: The local roughness exponent is found using the av- 
eraged wavelet coefficient (AWC) method. The solid line cor- 
responds to £zoc = 0.72, see Eq. J5J- This analysis was done 
on 237 lattices of size L — 256 with Daubechies wavelet of 
order 12. 

We apply the AWC method and the local window 
method to fracture profiles of lattice size L = 256. We 
measured a local fracture roughness of Q oc = 0.72 ± 0.02 
for the AWC method and Q oc = 0.70 ± 0.02 for the local 
window method, see Figs. |3 and [3J Both of these meth- 
ods underestimate the roughness exponent, with around 
0.03 for the AWC method and around 0.05 for the local 
window method. This was checked by creating artificial 
surfaces with known roughnesses using a wavelet method 
for generating the surfaces j2^|. For the local window 
method the deviations we found for the measured rough- 
ness exponent is consistent with the one found by Schmit- 
tbuhl et al. [2^. For the AWC method the deviation is 
consistent with the one found by Simonsen et al. [27j We 
therefore find the difference between £ and Q oc to be close 
to zero, which is a sign that anomalous scaling, as found 
by Zapperi et al. |3pJ in the 2D fuse model, is not found 
in the 2D central force model, and that there is only one 
roughness exponent, £ = 0.75 ±0.03, for the central force 
model. 

Anomalous scaling may be interpreted in the following 
way. Assume that the fracture process has produced a 
fracture of length 77 < L. The roughness measured over 
a window of size I < r\ will be 

w (l) = A(t]) , (6) 

where the prefactor A(r]) depends on 77 as a power law, 

A( V ) ~ , (7) 

so that w{l = rj) ~ 77^. Imagine we now fix a window size 
/ < 77 and follow the fracture roughness w(l) as 77 grows. 
The roughness within the window changes because the 
prefactor J7J) changes. This is not possible unless the 
cause of the anomalous scaling has its origin in statistical 
non-stationarity jsj, that is, the larger 77 is, the larger 
the sample over which the window I is applied to, and 
this will change the average on which the prefactor is 
build. Using the wavelet-based method proposed in Ref. 
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FIG. 3: Local fracture roughness measured by the local win- 
dow method, described in Eq. ||3J, where the averaged are 
done over the window size I for 237 lattices of system size 
L = 256. The solid line is Cloc = 0.70. 

[3l"| . we were unable to detect any such non stationarity, 
and as a consequence, ruling out anomalous scaling. 

When the disorder is broad, the correlation length for 
the percolation-like process £ is much greater then the 
system size L. Thus we can measure the correlation 
length exponent by using finite size scaling for the density 
of broken springs at fracture 

C 

Peff=Pc-jTjZ- (8) 
where p c is the critical density at which the network 



breaks down in the limit L — > oo, and p e ff is effective 
critical density at which the network breaks down for a 
finite L. This implies that the fluctuations in the density 
of broken springs at fracture will scale as 

"(Pe//) = ((Ptff) ~ (Peff} 2 ) 1/2 OC LT X I V . (9) 

Using a disorder of D = 20 which gives flat damage den- 
sity profiles, indicating that there is no localization in the 
fracture process, we obtain \jv = 0.57±0.10, see Fig. 01 
This corresponds to a value for v equal to 1.7±0.3, which 
by Eq. (JTJ gives C £ {0.73, 0.80}, which is consistent with 
the direct measurement of £. The value we found for v is 
different from the rigidity percolation value of 1.16 found 
by Moukarzel and Duxbury [22|. 

To conclude: We have studied the central force break- 
down model to test Eq. Q and found support for it mea- 
suring the roughness exponent of £ = 0.75 ± 0.03. The 
difference in global and local roughness exponent is small 
and a check for anomalous scaling of the wavelet coeffi- 
cients show no sign of such. The correlation length expo- 
nent for breakdown process was found to be v — 1.7±0.3 
giving consistent £ values using Eq. Q). The value for 
C and v are close to those of the fuse model, suggesting 
that these exponents might be similar for these models. 

This work was partially financed by the Norwegian Re- 
search Council and the computations were done under 
the auspices of the Norwegian HPC program, NOTUR. 
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